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Abstract 

We determine the position space fermion propagator in three dimensional QED based on Ward- 
identity and spectral representation. There is a new type of mass singularity which governs the long 
distance behaviour.lt leads the propagator vanish at large distancemore strongly than the mass term 
does. This term corresponds to Dynamical mass. Momentum space propagator is compared with the 
analysis of Schwinger-Dyson equation and our solution contains a non-perturbative effects beyond 
the quenched approximation with bare vertex. 
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I. INTRODUCTION 

In the previous paper we find that the propagator has a new type of mass singularity 
in QED3[l].The two kinds of self-energy are derived by Ward-Identity in evaluation of the 
second order spectral function. To include intermediate states with infinite numbers of pho- 
ton we exponentiate the one-photon matrix element. Usually three dimensional QED is 
thought as the high temperature limit of QED in the imaginary time formalism at finite 
temperature [2]. It is assumed that the three dimensional analysis leads the leading order 
of high temperature expansion results in four dimension [3]. In this respect first we study 
the position space structure of the propagator especially in the Yennie gauge. In this gauge 
we can avoid linear infrared divergences. In three dimension propagator has a logarithmic 
singularity even in the Yennie gauge. Fourier transformation to momentum space is done 
numerically. Low and high energy behaviour is consistent with other analysis in Euclidean 
space[6].In section II Bloch-Nordsieck approximation in three dimension is reviewed. Section 
III is for non-perturbative effects. In section IV mass singularity in four and three dimen- 
sion is compared and numerical analysis in momentum space is given. Main results are the 
followings. 

1 Based on Ward-Identity and dispersion-theoretic representation we find 
position space propagator Sp(x). 
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2 Gauge dependent terms are second-order corrected mass and linear infrared, 
divergences which vanishes in the Yennie gauge and others are gauge invariant. 

3 Second-order correction to the free propagator yields two kinds of gauge 
invariant mass singulairy;position dependent mass and self energy.The threshold 
behaviour is determined by self energy and dynamical mass is due to position 
dependent mass. Exponentiation of one-photon matrix element we obtain the full 
propagator. There exists dynamical symmetry breaking at e 2 /8iim = 1 for finite 
infrared cutoff. 



II. BLOCH-NORDSIECK APPROXIMATION TO QED 3 

First we consider about the charged praticle which emitt and absorb massless pho- 
tons. Usually this process was described by spectral function transition probability of par- 
ticle into particle and photon state. Multi-photon emitted from external line is introduced 
by ladder type diagrams which satisfy Ward-Identity. Let us beguin by dispersion theoretic 
description of the propagator [1] 



S F (p) = J d 3 xexp(ip-x) (n\TiP(x)ip{0)\n) 



j p z — in 1 + %t 

Sf(x) — J cIujSf(p,uj). (2) 

The field ip is renormalized and is taken to be a spinor with mass m .Here we introduce 
intermediate states that contribute the spectral function 

a(p 2 ) = (2it) 2 J2 53 (p-Pn) [ d 3 xexp(ip- x) (n\ij(x)\N) (N\^(0)\n) . (3) 



N 



Total three- momentum of the state \N) is p^.The only intermediates iV contain one spinor 
and an arbitrary number of photons. Setting 

\N) = \r;k 1 ,...,k n ), (4) 
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where r is the momentum of the spinor r 2 = m 2 ,and ki is the momentum of ith soft photo 
have 

<v> = / £ h. x '/ l^'^'* 2 ' D-'fr - r - E *o 

x (fi|^(x)|r; fci, .., fc n ) (r; fc ls ../c n |^(0)|fi) . 

Here the notation 

(/(*))o = 1, 

n 

(/(*))»= 

i=i 

has been introduced. Here we define matrix element 

T n = (n\i>\r;h,..k n ), 

T n=~ J d 3 xex P (tk n ■ x) (ft|7W|r; h, ..k n ^) , 

provided 

D X TMM = TijD x A,(x) = T^i-j^x) + d£(d ■ A{x))). 
T n satisfies Ward-Identity: 

d^T^x)) = -ei/>(x), 
d^j^x)) = e$(x), 
knnT^(r, k±, ..k n ) = eT n _i(r, k±, ..k n —±), r = m . 

Using LSZ the one photon matrix element is given 

7\ = (in\T(i/> in (x),ie J d 3 xlp in (y)^^ in (y)Af n (y))\r;k in^ 

= ie J fytzS* - V)l ^ ( y - ,) exp«* . „ + r ■ ,))✓<*, X)U(r, .) 

= 7 i A 2 5 -7M eM (^ A ) ex PW /c + r ) ■x)U(r,s), 

(r + k) A — m z 

where U (r, s) is a four-component free particle spinor with positive energy 

X^TTf \TT( \ l' r + m 

2^U{r, s)U(r, s) = ^ . 
s 
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The spectral function a is given by exponentiation of one-photon matrix element, 



T n = n T i(%)> 



(14) 



which yields an infinite laddar approximation 

— — exp(ir ■ x) exp(F), 



alx 



(15) 



(2tt) : 



exp(ifc • x)5(k 2 )9(k ) 



x e 2 tr 



(r + k) • 7 (r + k) • r u r • 7 + m 



(r + k) 2 — m? (r + k) 2 — m? 



2m n ^ 



_fk_ 

W) 2 



exp(ik ■ x)S(k 2 )9(k ) 



x e 



m 2 1 £(fc 

j A; 2 



2\1 



+ 



+ {d 



(16) 



(r • k) 2 r ■ k 

First we take the trace of T\T^~ for simplicity in the infrared. In this case we take the scalar 
part of them and assume pi(cu) = p2(cu).ln general case there are two kinds of spectral 
function. To aviod infared divergences we introduce photon mass p as an infrared cut off. It 
is helpful to use function D + (x) 

1 



D + (x) 



(2n) 2 i 
1 



exp(ik ■ x)d 3 k6(k°)5(k 2 - p 2 ) 

ixkdk exp(— px) 



J (kx) 



2^fW+ 



P A 



(2vr) 2 ^ 7 

to determine F.If we use parameter trick 

daexp(i(k + ie) • (x + ar)) = 
lim / adaexp(i(k + ie) ■ (x + ar)) = 



lim 



8nix 



exp(ik ■ x) 

k ■ r 
exp(ik ■ x) 
(k ■ r) 2 



(17) 



(18) 
(19) 



the function F is written in the following form 



F = ie 2 m 2 j adaD + (x + ar, p) — e 2 j daD + (x + ar, p) — ie 2 (d — l)—D + (x, p) 



e 2 m 2 exp(-px) 



87ir 2 



p 



+ x Ei(l, px)) + (d — 1)- — exp(— px), 

OTip 



(20) 
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where the function Ei(n, fix) is defined 

E iM = f ?*t^<l t . (21) 

It is understood that all terms which vanishes with // — > are ignored. The leading non 
trivial contributions to F are 

Ei(l, fix) = -7 - \n(nx) + O(fix), (22) 



Fl = ^(-±+z(l- ln(^) - 7 ) ) + O(fi), 



e 2 



F 2 = — (ln(/ix)+ 7 ) + 0(/i), 
F 9 = £(i-x)(rf-l) + 0(/x), (23) 

F = - 2) + £ + &r to <"*> - lnM - £ l(<i - 2 + ^ (24) 

where 7 is Euler's constant. Using integrals for intermediate state for on-shell fermion 

/ tttt? exp(ip • x) / d 3 r5(r 2 - m 2 ) = — — — -, (26) 
J (27r) J J m 2 + p z 

we set r = m in the phase space integral ; 

/d^ x f in 

- — —exp(ip-x) / =d 2 r exp(ir ■ x) exp(F(m, x)) 

[2tiY J v / r 2 + m 2 

o? 3 x ,. , exp(— mx) . . . . 

■exp(^p-x) exp(^ (to, x)). (27) 



/ 



(2tt) 3 y x 

Since we evaluated the matrix element F by the on shell limit of fermion and photon, 
F can be understood to describes self-energy [5]. In section IV and V we discuss dynamical 
mass,the renormalization constant ,and bare mass in connection of each terms. After angular 
integration we get the propagator 

a ( p ) = — dxsm(px)exp(-(m + B)x)(fix)- Cx+D , (28) 
2npJ 
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where 



1 (d-2) + ^-,B = ^(d-2 + 1 ),C=^,D=-^-. (29) 



87r/i 8iTm 87r 8tt 8imi 

Here we show up to 0(e 2 ) propagator a(p) 



m 



POO 

/ dx sin(px) exp(— mx) [A — Bx — Cx \n(fj,x) + D \n(fxx)} 
Jo 



2vrp J 
m 2 + p l [m A + p 1 ) 2 



where I\ , I2 are the following integrals 

h= r ^)-M-mx) 

Jo P 



-7 ln((m 2 + p 2 )/fi 2 ) ln((m - yf=p)/{m + ^f)) 



m? + p 2 2(m 2 + p 2 ) m 2 + p 2 

sin(pa;) exp(— mx) 



(31) 



p 

— m 



-x\n{jj l x)dx 



, 2 [ln((m- v ^)/(m + v ^)) + ln((m 2 +p 2 )/ /U 2 )-2(l- 7 )]. (32) 
[m 2 + p 2 ) 2 

Notice that the Yennie gauge (d = 2) is free from linear infrared divergences but associates 
a logarithmic divergences. In four dimension logarithmic divergence disappear and we have 
a free propagator in this gauge. This point is a clear difference between three and four di- 
mensional case. For defmitteness hereafter we choose the Yennie gauge. In this gauge the 
position space propagator is 

S F (x) = - f(i-yd + u f Xl> (~ UX K (u)du), (33) 
J Attx 

and it is written in momentum space 

S F (p) = / du — , = a(p) + trip), 

J p 2 — uj l + xe p 

a(p) — J d 3 x exp(ip ■ x)a(x), (34) 
a(x) = rneM-(jn + B)x M _ Cx+D ^ 

Finally we show the general spin dependent spectral function in the 0(e 2 ) for d = 2 gauge 

e 2 e 2 e 2 

F — (7 • r + ra)[— (—72; - x\n(mx) + (ln(mx) + 7)] - 7 • 7- — . (36) 

8n 871m 8nm 2 x 
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The first term is the same as the preavious one. The second term is independent of the cutoff 
/i and not significant in the infrared. In Euclidean space if we exponentiate it that becomes 

/mr 2 dr .. . exp(— mx) ir • 7 . ..ir-7 
exphr-x) — — -(1 + - _ )exp(( ; + 1)F) 

/mr 2 dr exp(— mrr) ir ■ x ,F , ,F a ■ r F 

TWT^ eMn ' ' X) 2vrx (1 + 2 )(C ° Sh( 2 } + ^ ( 2 )} - 

(37) 

This form corresponds to the infinite ladder graph of the fermion propagator with renormal- 
ized mass m and corrected vertex.In section IV we discuss physical meanings of each terms 
in a(x) and the dynamical effects in momentum space numerically. 

III. VANISHING BARE MASS AND VACUUM EXPECTATION VALUE OF 
CONDENSATE 

In this section we examine the renormalization constant and bare mass and study the 
condition of vanishing bare mass based on spectral representation. The spinor propagator in 
position space is expressed in the fllowing[3] 

S F ( X ) = - J{i 7 ■ d + cu) eW ^ X K (u)du (38) 

S F (p) = 1^a(p) + a(p). (39) 
V 

First term is related to wave function renormalization and second term is a mass func- 
tion which contains renormalized and dynamical mass. The equation for the renormalization 
constant in terms of the spectral functions read 



hm V(7-P + mo) = ^ r + (4Q) 
p^oo p 2 — m 2 + ie p^o J p 2 — uj 2 + ie 

Instead we determine them directly by taking the high energy limit of S F 



m Z^ = hm p 2 hr[S F (p)], (41) 

p— >oo 4 

Z^= hm ]tr[>y-pS F (p)]. (42) 



To determine Z 2 1 , first we show free case 



Sf W = (!7 ■ d) W^P = M^W = -7^^ (43) 
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where the dimension of the Sp\x) is equal to 1/x 2 . In momentum space 

] -x D ~ l dx ~ p~ 3 - D , (44) 



1 / ,00 a;sin(pa;)_ D _ lj _ 



V Jo V 
and this shows the ordinary expression 

1 



In this way we obtain 



lim -tr(7 • pS F (p)) ~ v ■ (45) 



Z 2 - ~ lim p-" = I (0<D) | 

\ 1 (0 = D) 



In the same way to evaluate Z 2 1 ,bare mass reads 

lim -p 2 tr(SF{p)) = m lim p 2 I 



xsin(pa;) D _, 
— -x ax 



px 

mp" D , (47) 



in the high energy limit if we use the formula 



m Z 2 l ~ m lim = , (48) 



and we obtain 

m ^-lZ^ > 0. (49) 



m 

Usually mass is a parameter which appears in the Lagrangean.For example chiral symmetry 
is defined for the bare quantity. In ref[4] the relation between bare mass and renormalized 
mass in the Schwinger-Dyson equation is discussed based on renormalization group equation 
and shown that the bare mass vanishes in the high energy limit even if we start from the finite 
bare mass in the theory. It suggests that symmetry properties must be discussed in terms 
of renormalized quanties.In QCD bare mass vanishes in the short distance by asymptotic 
freedom. And Dynamical mass vanishes too. In our apprximation this problem is understood 
that the two kinds of mass are always generated in which the second order terms in the 
coupling constans is gauge dependent but another is gauge invariant. In our position space 
propagator we can easily see the differences between bare mass and dynamical mass. But in 
momentum space its relation is complicated. In this case the symmetry of massless fermion 
is realized in the short distance. Condition of the vanishment of the bare mass depends on 
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the ratio of the renormalized mass and fixed coupling constant. There is a chiral symmetry 
at short distance where the bare or dynamical mass vanishes but its breaking must be 
discussed in terms of the values of the order parameter. Therefore it is interesting to study 
the possibility of pair condensation in our approximation. The vacuum expectation value of 
pair condensate is evaluated 



-trSp(x) = —2m lim 



exp(— mx) 



x 



(fix) 



-Cx+D 



( 



\ 



\ 



(50) 



(1 < D) 
finite (1 = D) 
oo (1 > D) J 

for fixed fi.We do not see confinement (Z 2 ^ 0) but condentation(finite (ipip)) occurs at 
D = Lit is understood that in QCD the anomalous dimension of quark or D in our model 
vanishes in the high energy limit thus Z 2 l =finite.In that case the vanishment of the bare 
mass is automatically satisfied[4].In the weak coupling limit we obtain Z 2 = l,m = m 
and (ipip) = oo . If we introduce chiral symmetry as global SU(n)i x SU(n) R ,it breaks 
dynamically into SU (n)y as in QCD [3] for D = l.Our model is also applicable to relativistic 
model of super fluidity in three dimension. 



IV. MASS SINGULARITY IN MOMENTUM SPACE 

In this section we study the effect of position depedent mass, self energy in momentum 
space 

M(x) = ^Hfxx), (51) 
e 2 

self energy = — In(fix). (52) 
The position space free propagator 

S F {x,m ) = -(ry • d + m ) (53) 

47TX 

is modified by these two terms which are related to dynamical mass and wave function 
renormalization.To see this let us think about the effect of spectral function in position 
space propagator 

exp(-(m + £ 
ai * ) = Anx ^ ] ^ M 8 ™- ( 54 ) 
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It is easy to see that the probability of particles which are separated with each other in 
the long distance is supressed by the factor (fxx)~ Cx ,and the self energy modifies the short 
distance behaviour from the bare 1/x to l/x 1 ~ D .The effect of self energy for the infrared 
behaviour of the free particle with mass m can be seen by its fourier transform [1] 

x2 sin(px)exp(-m.) = ^ + na + ^-./a, 

px Atcx 2 

~ (a/— p 2 — m) _1 ~ a near p 2 = m 2 . (55) 

Usually constant D represents the coefficent of the leading infrared divergence for fixed 
mass in four dimension. Therefore self energy has the same effects in three dimension as in 
four dimension [5] .In the non-linear Dyson-Schwinger equationof the fermion propagator or 
operator product expansion we can evaluate the dynamical mass which depends on momen- 
tum^] .Here we show that M(x) satisfies the same property as the dynamical mass. Let us 
define fourier transform of the scalar part of the propagator 

exp(-ms) , yCx+Ds = aE{p) = JIL f sin ( px ) exp(-mx)(fix)- Cx+D dx. (56) 
Airx 2%p J 

The momentum dependence of <7e(p) agrees well to the numerical solution of the Dyson- 
Schwinger equation with vanishing bare mass in the Landau gauge. It dumps as p~ 2 ~ D and 
stays constant at small p[6].It is easy to see that short distance behaviour is descrived by wave 
renormalization part (fix) D and the long distance behaviour is given by exp(— M(x)x) = 
(lix)~ Cx .The effect of (fj,x)~ Cx , first we check the o(e 2 ) contribution 

e 2 

0-M(p) = -^ J 2> ( 57 ) 



. 2 sin(pa;)exp(-mx) 

h = x xm((j,x)ax 

px Aixx 



o 



-m :W !5 1 v3 )+h( EM ) _ 2(1 _ 7)| , m 



Thus a secon order correction of the propagator dumps as (0(1) + ln(p 2 //x 2 ))me 2 /p 4 at high 
energy. It is not easy to see the effect of the factor (fix)~ Cx in momentum space. Therefore 
the scalar part of the propagator mj (m 2 +p 2 ) receives an additional correction by {jix)~ Cx .If 
the dynamical symmetry breaking and mass generation occur it has been discussed that the 
propagator will have a branch point on the real p 2 axis [7]. However the analysis in Minkowski 
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space is difficult due to the fact that the function (fix)~ Cx cannot be integrated analyti- 
cally. Therefore we see only the perturbative resutls.For this purpose it may be helpful to 
study the minkowski space by replacement from Euclidean distance R = Vx 2 — * iT in 

<j(x)[8], 

a(T) c expHmT) ( ^ r) _ tCT+D = exp(-imT) ^ct+o exp g j( _ {CT + D)) (5Q) 



exp(f(CT + Di))^T) D 



4tt 

provided 



cos(CT ln(/zT) + mT) - i sin(CT ln(//T) + mT) , (60) 



(ifiT)~ lCT+D = exp(-(C + Di))(fiT) D (cos(CT\n(fiT)) - i sm(CT \n(/iT))) . (61) 

Here we notice that cr(T) becomes complex valued function. Fourier transform into energy is 
given 

a(E) = lim 1- f cos(ET)a(T)dT, (62) 

= - J™ ^f* J T m°(T)dT. (63) 
In free case integration over space gives 



U J ArrV-t 2 + x 2 " 2;r/// 

1 2tt sin(mT) 



/ 



a (E)dE = lim 



t^oo 2T 7rm 2 
Thus energy distribution for free fermion becomes 



a Q (E) = ^ j ^cos(Et)a (t)dt 

m sin(mT) cos(ET) + £ cos(mT) sin(.ET) 



1ixTm{E 2 -m 2 + ie) 

Taking imaginary part 



(65) 



^ sin(2mT) 
9fro(#) = 2r ( 66 ) 

gives the residue at E 2 = m 2 .In free case imaginary part is 5(E 2 — m 2 ) and the continium 

contribution is given by an integral of the principal part in the region given by Q[E 2 — m 2 ).By 

numerical analysis using Maple we find that cr(E) osscillates and diverges in the region 

2 

E ~m eff = m + ^- + C (67) 

07T 
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V. SUMMARY 



Infrared behaviour of the propagator was examined in the Bolch-Nordsieck approximation 
to QED 3 . Even in the Yennie gauge there is a logarithmic infrared divergence but free 
from linear infrared divergences. Position dependent mass and self energy are both gauge 
independent in our approximation. Therefore in three dimension we examine the structure 
of the propagator in this gauge. In our approximation confinement property is given by 
position dependent mass. We evaluated the renormalization constant and bare mass based 
on spectral function sum rule and dynamical mass in momentum space is analised with 
finite infrared cutoff. In our approximation there seems to be a critical coupling constant 
for vacuum expectation value (ipip) which is due to wave function renormalization. This is a 
clear difference between lowest laddar approximation to the Schwnger-Dyson equation with 
bare vertex and ours. In Euclidean space we see the momentum dependence of the propagator 
which has been known in the anlysis of the Schwinger- Dyson equations except for the wave 
function renormalization [6]. 
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